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ABSTRACT 

I apply the well controlled Hydro-PM approximation of Gnedin & Hui to model the 
column density distribution of the Lyman-alpha forest for 25 different flat cosmolog¬ 
ical scenarios, including variants of the standard CDM, tilted CDM, CDM with a 
cosmological constant, and CHDM models. I show that within the accuracy of the 
HPM approximation the slope of the column density distribution reflects the degree 
of nonlinearity of the cosmic gas distribution and is a function of the rms linear den¬ 
sity fluctuation at the characteristic filtering scale only. The amplitude of the column 
density distribution, expressed as the value for the ionizing intensity, is derived as a 
function of the cosmological parameters (to about 40% accuracy). The observational 
data are currently consistent with the value for the ionizing intensity being constant 
in the redshift interval z ~ 2 — 4. 

Key words: cosmology: theory — intergalactic medium - quasars: absorption lines 
- methods: numerical - hydrodynamics 


1 INTRODUCTION 

Recent years witnessed a rapid development in our under¬ 
standing of the Lyman-alpha forest, a phenomenon, poten¬ 
tially offering cosmologists the most accurate probe of the 
Intergalactic Medium (IGM) at intermediate cosmological 
redshifts (z ~ 2 — 4). On one hand, observations with the 
Keck telescopes (see, for example, Hu et al. 1995; Lu et al. 
1996; Cristiani et al. 1996; Kirkman & Tytler 1997; Kim et 
al. 1997; and D’Odorico et al. 1997) provided an impressive 
amount of high quality data. On the other hand, cosmo- 


logical hydrodynamic simulations (Cen et al. 1994 

; [Zhang, 

Anninos, & Norman 19951; 

Hernciuist et al. 1996; 

Miralda- 

Escude et al. 1996 

; Wadsley & Bond 1997 

Zhang et al. 1997) 


offered a theory whose a priori predictions can now not only 
be quantitatively compared to the observational data, but 
also, and perhaps even more importantly, be used to con¬ 
straint the cosmology with the Lyman-alpha forest observa¬ 
tions. 


While a possibility still remains that the true nature 
of the Lyman-alpha is still uncovered, or that there exist 


a genuine population of individual absorbers (Bahcall & 

Salpeter 1965|; Arons 1972; Black 1981; 

Ostriker & Ikeucli: 

1985; 

[keuchi & Ostriker 19861; Rees 1986; Ikeuchi 1986; |Reef 

1988; 

Bond, Szalay, & Silk 1988), this paper is solely based 


on the theory that Lyman-alpha absorption arises from the 
fluctuating IGM, and the fluctuations in the IGM are small- 
scale siblings of larger-scale structures, that one sees in the 


© 0000 


distribution of galaxies in the sky, as argued for by cosmo¬ 
logical hydrodynamic simulations. 

However, cosmological hydrodynamic simulations may 
be the worst theoretical tool one can imagine. They re¬ 
quire enormous (compared to the vast majority of theoreti¬ 
cal methods) computational resources, and at the same time 
demand high expertise from a researcher, who ought to un¬ 
derstand a vast number of subtle numerical and modeling 
biases arising from using such a complicated tool as a cos¬ 
mological hydrodynamic code. The development of the code 
takes months, if not years, of human labor, and currently 
there exist only a few groups who engage in modeling the 
Lyman-alpha forest with hydrodynamic simulations. This 
situation will change in the future, when supercomputers are 
faster and more widely available, and when standard hydro- 
dynamic codes become publicly accessible, but the current 
situation vows for the use of approximate methods, which 
sacrifice the accuracy for the sake of simplicity and speed. 
The simplicity and speed are especially important because 
currently there exist a large variety of cosmological models 
that need to be tested against the observational data; while 
one or two specific cosmological models can be best studied 
with the hydrodynamic simulations, it is currently beyond 
reason to investigate dozens of models with the hydrody¬ 
namic simulations. 


Up to date, two different approximations have been used 
for the purpose of model ing the Lyman-alpha forest: the log- 
normal approximation (Bi, Borner fc Chu 1992: 


Gnedin & 
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Hui 199 6; [tii fc Davidsen 1997 ) and the Z crdovich approxi¬ 
mation ( Doroshkevkh^&^Shandarin^977 ; McGill 1990|; Hui 


Gnedin, & Zhang 1997). While both approximations are very 
efficient and much easier to implement than a cosmological 
hydrodynamic simulation, they are uncontrolled approxima¬ 
tions, i.e. the range of their applicability and their accuracy 
is not well established. As a result, their predictions cannot 
be directly used in a quantitative comparison with the ob¬ 
servational data, since it is not clear what the error is to 
assign to the theoretical prediction based on an approxima¬ 
tion (since an approximation is used, it has to have an error 
associated with it). 

In this paper I use a different a pprox imation, called 
HPM, and developed in Gnedin & Hui (1998). In its essence, 


the HPM method uses a simple Particle-Mesh (PM) solver 
modified to account for the effect of gas pressure. Be¬ 
cause there exist a tight correlation between gas temper¬ 
ature/pressure and gas density (“equation of state”) in the 
low density IGM in the mildly nonlinear regime (cosmic 
overdensity 8 S-) 10), it is possible to compute the gas tem¬ 
perature and pressure at every point directly from the value 
of the cosmic gas density at this point. Thus, there is no 
need to introduce a special equation for the gas tempera¬ 
ture as in a full hydrodynamic solver. As a result, the HPM 
approximation is only about 25% slower (due to the over¬ 
heard of computing the equation of state) than a simple PM 
solver, and substantially faster than a full hydrodynamic 
solver (due to both fewer computations at each time-step 
and fewer time-steps), while delivering results which are ac¬ 
curate to about 15% in the point-by-point comparison with 
a full hydrodynamic simulation of precisely the same cosmo¬ 
logical model (for 8 Sa 10). 

Standard (i.e. unmodified to account for the gas pres¬ 
sure) PM simulations have been also used to address prop¬ 


erti es of the Lyman-alpha forest (Petitjean, Miicket fc Kate; 


1995; Miicket et al. 1996; Weinberg 1997). While PM sim¬ 
ulations are most likely more accurate than lognormal or 
Zel’dovich approximations, they still cannot quite reach the 
level of accuracy achieved by the HPM metho d. For exam¬ 
ple, as has been shown by Gnedin & Hui (1998), a standard 
PM simulation predicts about a factor of 3 lower densities 
than the full hydrodynamic simulation; it is possible to im¬ 
prove upon the PM simulation by appropriate smoothing of 
the initial conditions, but even in this case a PM simulation 
is only able to reproduce the gas density to about 30-40% 
for 8 ^ 10, whereas HPM reaches a 15% accuracy. 

In general, it does not seem to be much rationale in us¬ 
ing the standard PM simulation instead of HPM, because 
25% increase in the CPU time in the HPM approach com¬ 
pared to simple PM buys at least a factor of 2 in the error 
reduction. Also, from a technical point of view, it does not 
require much human effort to convert a PM code into an 
HPM code. 

Having chosen the tool to study the Lyman-alpha for¬ 
est, it is most natural to start with investigating the column 
density distribution, since it is the quantity best known ob- 
servationally (to within about 10-15%). Since the main ad¬ 
vantage of using the HPM approximation is its speed, I con¬ 
centrate in this paper on analysing the properties of the 
column density distribution of the Lyman-alpha forest for 
a large set of different cosmological models (25 altogether) 
with the main goal to understand what physical parameters 


the column density distribution depends upon. Testing of 
specific cosmological models is only a secondary goal of this 
paper. 

I begin my analysis with discussing the models and 
characteristic scales in Section 2. In the same section I also 
discuss another approximation I have to adopt to conclude 
my study, and carefully investigate the total error induced 
by using approximations rather than an exact method (in 
this case a cosmological hydrodynamic simulation). The re¬ 
sults are presented in Section 3, and Section 4 concludes the 
paper with the discussion. 


2 METHOD 


2.1 Cosmological models 


With the help of HPM, it is possible to simulate a large 
set of cosmological models with only modest computational 
resources. Table [I] lists all 25 simulations considered in this 
paper. Simulations have been performed with 256 3 particles 
on a 256 3 mesh. An average 256 3 simulation takes about 40 
CPLT-hours on a R10K processor. In addition, 4 simulations 
have been rerun with 512 3 particles on a 512 3 mesh, to verify 
numerical convergence. Those four models are shown in bold 
face in Table |I|. Each of 512 3 simulations requires about 450 
CPU-hours on a R10K processor, and it is impractical at the 
moment to perform all simulations listed in Table [I] with the 
512 3 mesh. I therefore use the 512 3 simulations to investigate 
numerical convergence and accuracy of the results obtained 
with 256 3 models. 

Placement of a simulation box has been chosen so as 
to insure that the characteristic scale over which linear 
baryon fluctuations a re smoothed, the so called filtering scale 
( |Gnedin fc Hui 1998| ), is resolved with at least 3 cells. This 
dictates the comoving cell size of 10h -1 kpc for Ho = 1 mod¬ 
els and 15h -1 kpc for LCDM models, with the exception of 
the SCDM.2K and SCDM.2L models, whose low reionization 
temperature and late reionization (and therefore a small fil¬ 
tering scale) required to use the cell size of 7h -1 kpc. 

Cosmological parameters and labels for all models are 
summarized in Table [I]. Here Ho is the total matter density 
parameter, fif, is the baryon density parameter, Ha is the 
cosmological constant density parameter, and Q, v is the den¬ 
sity parameter in massive neutrinos. As usual, h denotes the 
Hubble constant in units of 100 km/ s/ Mpc, n is the power- 
law index of the primordial spectrum of density fluctuations 
(n = 1 is the scale-free Harrison-Zel’dovich spectrum), and 
os is the rms top-hat density fluctuation on 8h -1 Mpc at 
2 = 0 , 


oj = 


b 2 rtk 

— r P(k)Wi H {k x 8 h- 1 Mpc). 


In order to simulate the evolution of the IGM un¬ 
der the HPM approximation, an equation of state (i.e. the 
temperature-density relation) as a function of time ought 
to be specified. In general, the equation of state depends 
on the time evolution and spectral shape of the ionizing 
background, which is a two-dimensional function. In order 
to limit the possible choice of time-dependent equations of 
state to a manageable set, I adopt the power- law e quation 
of state in the form derived in Hui & Gnedin (1998), 
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Table 1. Cosmological Models 


Model 

Ho 

Hf, 

Ha 

H„ 

h 

n 

C8 

&34 

Ti, rei 

-2-rei 

7 

t 4 

a F 

°F,BOX 

COBE 

SCDM.1A 

1 

0.05 

0 

0 

0.5 

1 

0.8 

2.22 

2.5 

7 

1.49 

0.91 

2.39 

2.15 


SCDM.2A 

1 

0.05 

0 

0 

0.5 

1 

0.7 

1.95 

2.5 

7 

1.49 

0.91 

2.09 

1.88/2.00 


SCDM.2D a 

1 

0.05 

0 

0 

0.5 

1 

0.7 

1.95 

2.5 

7 

1.49 

0.91 

2.09 

1.88/2.00 


SCDM.2E a 

1 

0.05 

0 

0 

0.5 

1 

0.7 

1.95 

2.5 

7 

1.49 

0.91 

2.09 

1.88 


SCDM.2G 

1 

0.05 

0 

0 

0.5 

1 

0.7 

1.95 

1.5 

7 

1.51 

0.61 

2.23 

2.04 


SCDM.2H 

1 

0.05 

0 

0 

0.5 

1 

0.7 

1.95 

4.0 

7 

1.46 

1.32 

1.96 

1.74 


SCDM.2J 

1 

0.05 

0 

0 

0.5 

1 

0.7 

1.95 

2.5 

15 

1.61 

0.74 

1.92 

1.70 


SCDM.2L 

1 

0.05 

0 

0 

0.5 

1 

0.7 

1.95 

2.5 

5 

1.33 

1.21 

2.26 

2.00 


SCDM.3A 

1 

0.05 

0 

0 

0.5 

1 

0.6 

1.67 

2.5 

7 

1.49 

0.91 

1.79 

1.62 


SCDM.4A 

1 

0.05 

0 

0 

0.5 

1 

0.5 

1.39 

2.5 

7 

1.49 

0.91 

1.49 

1.35/1.43 


SCDM.5A 

1 

0.05 

0 

0 

0.5 

1 

1.2 

3.34 

2.5 

7 

1.49 

0.91 

3.58 

3.23/3.52 

V 

TCDM.2A 

1 

0.05 

0 

0 

0.5 

0.9 

0.83 

1.95 

2.5 

7 

1.49 

0.91 

2.09 

1.84 

V 

TCDM.4A 

1 

0.05 

0 

0 

0.5 

0.7 

1.16 

1.95 

2.5 

7 

1.49 

0.91 

2.10 

1.74 


LCDM.1A 

0.35 

0.05 

0.65 

0 

0.7 

1 

0.94 

1.85 

2.5 

7 

1.53 

1.25 

1.94 

1.69 

V 

LCDM.1C 

0.35 

0.05 

0.65 

0 

0.7 

1 

0.94 

1.85 

1.5 

5 

1.43 

0.96 

2.26 

2.06 

V 

LCDM.2A 

0.35 

0.05 

0.65 

0 

0.7 

0.96 

0.80 

1.49 

1.5 

5 

1.43 

0.96 

1.80 

1.62 

V 

LCDM.3A 

0.35 

0.03 

0.65 

0 

0.7 

1 

1.04 

2.14 

2.5 

7 

1.51 

1.02 

2.27 

2.01 

V 

LCDM.4A 

0.35 

0.03 

0.65 

0 

0.7 

0.95 

0.85 

1.62 

2.5 

7 

1.51 

1.02 

1.72 

1.50 

V 

LCDM.5A 

0.4 

0.036 

0.6 

0 

0.65 

1 

1.02 

2.12 

2.5 

7 

1.51 

1.04 

2.25 

2.00 

V 

CHDM.1A 

1 

0.05 

0 

0.1 

0.5 

1 

0.91 

2.13 

2.5 

7 

1.49 

0.91 

2.26 

2.02 

V 

CHDM.2A 

1 

0.05 

0 

0.15 

0.5 

1 

0.85 

1.69 

2.5 

7 

1.49 

0.91 

1.78 

1.57 

V 

CHDM.3A 

1 

0.05 

0 

0.2 

0.5 

1 

0.81 

1.35 

2.5 

7 

1.49 

0.91 

1.41 

1.22 

V 

CHDM.4A 

1 

0.05 

0 

0.1 

0.65 

1 

1.21 

3.29 

2.5 

7 

1.50 

0.99 

3.46 

3.12 

V 

CHDM.5A 

1 

0.05 

0 

0.2 

0.65 

1 

1.11 

2.10 

2.5 

7 

1.50 

0.99 

2.19 

1.89 

V 

CHDM.6A 

1 

0.07 

0 

0.1 

0.6 

0.9 

0.72 

1.55 

2.5 

7 

1.52 

1.09 

1.61 

1.41 

V 


a A different random realization of the SCDM.2A model. 


T(p) = T 0 (l + 8) 


7-1 


(i) 


where 5 is the cosmic overdensity, and To and 7 are func¬ 
tions of redshift. Speci fically , I adopt the evolution of To and 
7 from Hui & Gnedin (1998), equations (17) and (19). In this 
case the evolution of the equation of state is determined by 
only two parameters, the epoch of sudden reionization z re i, 
and the temperature of the gas just after reionization, T re i. 
Those two parameters are shown in Table [I] for all mod¬ 
els, with the reionization temperature expressed in units of 
10 4 K, 


T4 tre i = T r ei/10 4 K, 


with reasonable values ranging from Tj_ re i ~ 1 to T 4 . re i ~ 4. 
For the purpose of this paper it is only important that equa¬ 
tion (i,lj|) gives an equation of state, and by changing param¬ 
eters T re i and z re i I can span a reasonable range in possible 
equations of state. As will be shown below, the specific form 
of the equation of state plays a relatively minor role in deter¬ 
mining the column density distribution of the Lyman-alpha 
forest. 

Given the specific form of the evolution of the equation 
of state, it is then possible to compute the filtering scale 
as a function of redshift for a given cosmological model. In 
general, the wavenumber k_p corresponding to the filtering 
scal e is given by the following expression (Gnedin & Hui 

199C] 


1 


k 2 F (t) 


D+(t) 







( 2 ) 


where D+(t) is the growing mode of linear fluctuations, and 


kj = — ^AnGg (3) 

cs 

is the wavenumber corresponding to the Jeans scale. Here 
cs is the sound speed of the cosmic gas, and g is the total 
cosmic mean density. For the fully ionized cosmic gas with 
the mean molecular weight g = 0.6 and with the equation 
of state given by (jl|), kj can be expressed as 

kj(z) = 8.44fl 1 ^ 2 /tMpc ~ 1 ( L5 ^ Z) ) ' . (4) 


where T 4 = To/10 4 K. In general, the filtering scale is a 
complicated function of time, but for a large set of reason¬ 
able equations of state it is about a half of the Jeans scale 
at 2 ~ 3, Uf ~ 2kj. It is therefore useful to introduce a 
measure of nonlinearity of the gas distribution of a given 
cosmological model at z = 3, a 34, defined as 

k 2 rlk 

-^rP(k, z = 3) exp(-2fc 2 /fc 3 2 4 ), (5) 

where 



fc 3 4 = 34STq4 2 Mpc 1 


( 6 ) 


is the wavenumber corresponding to one half of the Jeans 
scale for the cosmic gas with T 4 = 1 and 7=1. 5 at 2 = 3. 
The partic ular choice of filtering was discussed in Gnedin & 
Hui (|l99lj). The quantity 074 should therefore indicate an 
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approximate degree of nonlinearity of the gas distribution 
in a given cosmological model at z ~ 3. A more accurate 
indicator is of course the rms linear density fluctuation at 
the filtering scale, 


°f0) = 


k 2 dk 

2n 2 


P(k,z)ex p [— 2k 2 /k%(z)^ 


(7) 


and kp(z) is a function of time. While of is a more accurate 
measure, it is also more difficult to compute (since it requires 
evaluation of Icf). On the other hand, 034 is less accurate 
but easier to compute, and it can be useful for quick, order- 
of-magnitude estimates. 

Since a simulation is performed on a finite box, only a 
finite range of scales is present in the simulation. Thus, the 
rms density fluctuation computed in a simulation will be 
smaller than the one obtained using the full range of scales. 
To be able to account for this effect, I introduce another 
quantity, of,box, similar to equation (|7j), but with integra¬ 
tion performed over only a finite range of scales from k m i n 
to fc max which are present in the simulation box: 


o>, box (2) = 


k 2 dk 

2n 2 


P(k,z)ex p [—2 k 2 /kp(z)^ 


( 8 ) 


The small scale cut-off, fc max , is not that important because 
the exponential filtering insures that, as long as the fil¬ 
tering scale 1//cf is resolved with at least three cell sizes 
(/cf < 1/(31) where l is the cell size), any contribution to 
the integral at k ~ fc max = n/l (the Nyquist frequency) is 
suppressed by a factor exp[— 2 ( 37 r) 2 ], which is equal zero for 
any practical purpose. 

Values of 0 - 34 , ctf(z = 2.85), and o\f,box (2 = 2.85) are 
shown in Table ^ (for four models for which both 256 3 and 
512 3 simulations have been run, both values of of,box (2 = 
2.85) are shown). Also, for reference purposes, I show values 
for T 4 and 7 which define the equation of state of the IGM 
at 2 = 2.85. The redshift value of 2 = 2.85 corresponds 
to the epoch where most abundant observation data on the 


column density distribution exist (Hu et al. 1995 
& Tytl 


sr 1997 


Kirkman 


Kim_et al. 1997). The check mark in the 


last column of Table [l] indicates whether a particular model 
is CO BE normalized. The CO BE normalization was done 


using Bunn & White (1997) fits and assuming the gravity 


wave contribution for n < 1 models. 

I 11 order to specify initial conditions for a simulation, a 
linear matter transfer function needs to be computed. For 
my LCDM simulations I used transfer functions computed 
using the linear gravity code (existing in the literature fits 
are not sufficiently accurate due to the higher baryonic frac¬ 
tion in LC DM m odels compared to S2o = 1 models); Hu & 
Sugiyama ( 199C| ) approximation was used to compute trans¬ 
fer functions for SCDM and TCDM simulations, and Ma 


(1996) fits were used for CHDM simulations. 

Finally, a few words about the choice of models. Ap¬ 
parently, only COBE normalized models can be considered 
viable. The set of SCDM models, which are not COBE nor¬ 
malized, is therefore used as a kind of “training set”, with 
which dependence on different parameters is analysed in or¬ 
der to develop a general understanding of how the column 
density distribution of the Lyman-alpha forest depends on 
the underlying cosmology. 


2.2 The column density distribution 

The purpose of this paper is to investigate the column den¬ 
sity distribution of the Lyman-alpha forest, and I therefore 
need a way to compute the column density distribution from 
the output of a simulation. The standard way of accomplish¬ 
ing this is to generate sy nthetic spectra and then fit Voigt 
profiles (Dave et al. 1997). However, the Voigt profile fitting 


is extremely CPU intensive, and full Voigt profile analysis of 
results of one simulation would require more computer time 
than the simulation itself consumed. It is therefore feasible 
to use an approximate technique and then properly account 
for the error induced. 

Fortunately, there is such an approximate method called 
t he “Density-Peak An s atz” (hereafter DPA), introdu ced in 
( |Cnedin fe Hui 1996 ; Hui, Gnedin, & Zhang 1997). The 
DPA approximation is based on the assumption that a single 
Lyman-alpha line arises from a peak in the cosmic gas den¬ 
sity along the line of sight to a distant quasar. Then, given 
the value of the cosmic density p = 1 + 8 at the peak, and 
the second derivative of the density along the line of sight 
(the first derivative vanishes at the peak), a value for the 
column density can be associated with the peak as follows: 


jYh, = 1.45 x 10 13 cm" 2 T 4 “°- 7 


n b h 2 V/0.5 
0.0125/ VJ 21 


1 + 2 ' 


(1 + 5 ) 


2 - 0 . 7 ( 7 - 1 ) 


-(1-0.35(7-1)) 


d 2 ln(l + 5) 
dx 2 


- 1/2 


, 0 ) 


and x being measured in Mpc. Here J 21 is the ionizing in¬ 
tensity, 




V dv/v 


f Ojjdv/v 10 -21 ergcm -2 s -1 Hz 1 sr -1 ’ 

where J u is the radiation intensity and a v is the hydrogen 
photoionization cross-section. 

I 11 order to test the DPA approximation, I have com¬ 
puted synthetic spectra for random lines of sight through 
the simulation box for several models at 2 = 2.85 and fitted 
Voigt profiles to Lyman-alpha fore st lines using AU TOVP 
automatic Voigt profile fitting code (Dave et al. 1997), kindly 
provided to me by Romeel Dave. In order to do that, I need 
to assume a value for the ionizing intensity J 21 , and also 
specify the signal-to-noise ratio S/N for Voigt profile fitting 
(the DPA column densities are independent of the adopted 
signal-to-noise ratio). Figure |l] shows the scatter plot of 
column densities from the DPA approximation versus col¬ 
umn densities computed using Voigt profile fitting for the 
SCDM.2A model assuming J 21 = 0.17 (which gives a rea¬ 
sonable fit to the observational data from Hu et al. 1995) 
and S/N = 50. Note, that while there are absorption lines 
on which the DPA approximation fails severely, in large the 
agreement is remarkable, keeping in mind the simplicity of 
the DPA method. 

Column density distributions from the DPA method 
and Voigt profile fitting for several models are shown in Fig¬ 
ure ^ with thin and bold lines respectively. The solid line 
shows the SCDM.2A model with J 21 = 0.17 and S/N = 50, 
while the dotted line shows the same model but with S/N = 
100. I also show the CHDM.3A and LCDM.3A models on 
the same plot, as well as the SCDM.2A model at 2 = 3.70. 

One can see that DPA agrees very well (to about 3%) 
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Table 2. Error budget 


N h ,(VP) 

Figure 1. Scatter plot of column densities computed using the 
DPA approximation vs column densities obtained by Voigt profile 
fitting for the SCDM.2A model with J 21 = 0.17. 



N hi (cm- 2 ) 

Figure 2. Column density distributions computed using the 
DPA method (thin lines) and Voigt profile fitting (bold lines) for 
the SCDM.2A model with J 21 = 0.17 and S/N = 50 (solid lines), 
the SCDM.2A model with J 21 = 0.17 and S/N = 100 (dotted 
lines), the CHDM.3A model with J 21 = 0.25 and S/N = 50 
(short-dashed lines), the LCDM.3A model with J 21 = 0.3 and 
S/N = 50 (long-dashed lines), and the SCDM.2A model with 
J 21 = 0.25 and S/N = 50 at z = 3.70 (dot-dashed lines). The 
last three models are plotted with the successive vertical offsets 
of 0.5 dex for clarity. 


Source of error 


Slope Amplitude 


HPM is not exact 3% 3% 

DPA is not exact 3% 2% 

Finite box size 13(3)% 18(5)% 


Total error 


14(5)% 19(6)% 


with the column density distribution from Voigt profile fit¬ 
ting for column densities between about 10 13 cm- 2 and 
10 14 cm -2 independently of the signal-to-noise ratio in the 
Voigt profile fitting procedure, but fails at lower and higher 
column densities. At the lower end, the DPA column den¬ 
sity distribution falls off due to the finite resolution of the 
simulations; in this regime the DPA approximation becomes 
unreliable, and, at the same time, the Voigt profile fitting 
technique becomes sensitive to the assumed signal-to-noise 
ratio in the synthetic spectrum. One can expect, however, 
that higher resolution (larger mesh size) simulations would 
improve agreement between the DPA and Voigt profile fit¬ 
ting for lower column densities. 

For column densities in excess of 10 14 cm" 2 , i.e. for satu¬ 
rated lines, the column density distribution computed using 
Voigt profile fitting shows a feature, while the DPA column 
density distribution continues as a smooth function of the 
column density. It is therefore plausible to presume that the 
feature in the column density distribution obtained by Voigt 
profile fitting is an artifact of the procedure rather than a 
change in intrinsic properties of absorbers. However, since all 
available observational data were analysed using the Voigt 
profile fitting technique, the same artifact is present in data 
as well. Thus, it seems reasonable to restrict the range of ap¬ 
plicability of the DPA method to column densities between 
about 10 12 ' 6 cm- 2 and 10 14 '°cm- 2 , where the agreement is 
excellent. As a side note I point out here that this range of 
column densities is also within the range of applicability of 
the HPM approximation, which gives accurate results only 
for overdensities S & 10. 


2.3 The error budget 

The HPM approximation for the hydrodynamics of the low 
density IGM, combined with the DPA approach, is a pow¬ 
erful and accurate method of computing the column density 
distribution of unsaturated Lyman-alpha absorption lines. 
Both approximation, the HPM, and the DPA, are controlled 
approximations, i.e. their range of applicability and their ac¬ 
curacy is well known. It is therefore possible to compute the 
accuracy with which the column density distribution can be 
computed using the HPM-DPA combination. 

It is now customary to approximate the column density 
distribution of the Lyman-alpha forest with the power law 
over a finite range of column densities. Thus, I can convert 
the error in the computed column density distribution to er¬ 
rors in the amplitude and the slope of the power law. There 
are three different sources of error in the HPM-DPA cal¬ 
culation, which are summarized in Table 0. First, the HPM 


itself is not exact. Analysis of Fig. 11 of Gnedin & Hui (1995) 
shows that the HPM introduces about 3% error to both the 
slope and the amplitude of the power law. Analysis of Fig. 
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Figure 3. Column density distributions for three variants of the 
SCDM model: erg = 0.5 (dotted lines), erg = 0.7 (solid lines), and 
erg = 1.2 (dashed lines). Thin lines show 256 3 simulations and 
bold lines show 512 3 simulations. For ag = 0.7 model three dif¬ 
ferent random realization with 256 3 resolution and two different 
random realizations with 512 3 resolution are shown. The square 
bracket shows the region where numerical convergence with re¬ 
spect to the box size is achieved. Lines of different type are offset 
vertically for the sake of clarity. 


also shows that a 3% error in the slope and a 2% error 
in the amplitude is introduced by using the DPA instead of 
Voigt profile fitting for column densities between 10 13 cm -2 
and 10 14 cm -2 . Finally, due to the finite size of a compu¬ 
tational box, only a subset of all possible initial conditions 
can be simulated (only one “realization”), and this statis¬ 
tical undersampling also introduces an error. In order to 
estimate the error due to the realization dependence, I have 
performed three different realizations of the same cosmolog¬ 
ical model with the 256 3 mesh and two realizations with the 
512 3 mesh (models SCDM.2A, SCDM.2D, and SCDM.2E). 

Figure ra shows column density distributions for all re¬ 
alizations, normalized as explained in §2.4. The difference in 
slopes and in power law amplitudes for three 256 3 realiza¬ 
tions reaches 10% and 18% respectively, while 512 3 simula¬ 
tions have an error a factor of three less in both the slope 
and the amplitude. Thus, performing 512 3 simulations for 
all 25 model would substantially reduce the theoretical er¬ 
ror, but, as I mention above, since this is beyond the prac¬ 
tically possible at the moment, the total error from Table || 
should be properly included in all comparisons between the 
simulations and the observations. 

While Table ^ gives a magnitude of the error, it does 
not determine the statistical property of this error. In other 
words, it does not specify what the probability is that a 
given model would produce results accurate to within the 
error from Table While there is no simple way to find out 
the probability distribution of the error, the error in Table ^ 
should be considered as about 95% confidence interval from 


Table 3. Column density limits 


Redshift 

1 §(Vh i,min ) 

lg(V H i,MAx) 

2.31 

12.6 

13.2 

2.85 

12.6 

13.6 

3.35 

13.0 

14.0 

3.70 

13.4 

14.0 


the way it was computed. This assertion is also supported by 
comparison between the observations and the simulations in 
the following section. 


2.4 Observational data 


I now turn to directly comparing the results of the simula¬ 
tions with the observational data. Recent Keck observations 


of the Lyman-alpha forest (Hu et al. 1995 Lu et al. 199f; 
Kirkman & Tytler 1997; Kim et al. 1997) give the column 
density distribution of the Lyman-alpha forest at four dif¬ 
ferent redshifts: z = 2.31, z = 2.8 5, z = 3.35, and z = 3.70. 
The data from Kirkman & Tytler (1997) have a median red- 
shift of z = 2.70 and about 10% higher amplitude than the 
data from Hu et al. ( 1995| ). Some part of this difference is 
due to the slightly different redshift of Kirkman & Tytler 
observations, and some other part (not quite certain at the 
moment) is due to the Cosmic Variance. But since the cur¬ 
rent error of my simulations is about 20% in the amplitude, 
both data samples are consistent with each other for my 
purpose. 

As Fig. [l]also demonstrates, only within a limited range 
of column densities there is a consistency between 256 3 and 
512 3 simulations, i.e. only with this final range a numeri¬ 
cal convergence is achieved. Thus, only this range of column 
densities can be used in comparing observations and simu¬ 
lations. In addition, this range is also a function of redshift, 
and a figure similar to Fig. at each value of redshift is 
used to determine the range of column densities at each 
redshift within which a numerical convergence is achieved. 
Whenever the range of convergence exceeds the range of ap¬ 
plicability of the HPM-DPA approximation, a smaller of two 
ranges is used (thus maintaining both the accuracy of the 
HPM-DPA approximation and the numerical convergence of 
simulations). As can be seen from Fig. the range of con¬ 
vergence is a weak function of the amplitude of the density 
fluctuations, and thus the same range of column densities 
can be used for all models. 

The specific limits are summarized in Table I only 
point out here that 512 3 simulations have not been contin¬ 
ued beyond z = 2.85, and the fitting interval at z = 2.31 
is determined by extrapolation of fitting intervals at higher 
redshifts. One thus has to keep in mind that numerical con¬ 
vergence has not been firmly established for simulation re¬ 
sults at z = 2.31. Therefore, I will only use simulation results 
at z = 2.85, z = 3.35, and z = 3.70 for quantitative com¬ 
parison with observations, and will use model predictions at 
z = 2.31 for illustrative purposes only. 

For a given range of column densities, I can fit a power 
law (simultaneously in the slope and amplitude) to both 
the results of simulations and the observational data. I then 
fix the amplitude of the column density distribution from 
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Figure 4. The slope of the column density distribution as a 
function of the rms linear density fluctuation at the filtering scale 
c tjp for all models at four different redshifts: 2 : = 2.31 (crosses), 
2 = 2.85 (circles), 2 = 3.35 (squares), 2 = 3.70 (triangles). Open 
symbols mark 25 256 3 simulations, and filled symbols show 4 512 3 
simulations. The solid line shows the fit [l^, and dotted lines show 
d= 13% deviation from the fit. 


the simulation by requiring that it coincides with the am¬ 
plitude fitted to the observational data at the specific value 
of the column density A^hi,fit = 3 x 10 13 cm 2 (for all four 
values of redshift). Thus, for a given simulation this proce¬ 
dure produces two numbers: the slope of the column den¬ 
sity distribution (3 and the amplitude of the column density 
distribution, which can be expressed as J 21 . It is the time 
evolution and the dependence on cosmological parameter of 
those two quantities that I concentrate on in this paper. 


3 RESULTS 


3.1 The slope of the column density distribution 


I now focus on investigating the dependence of the slope of 
the column density distribution (3 on cosmological parame¬ 
ters. 

Figure [I] presents the major result of this paper, the 
dependence of the slope of the column density distribution 
on the rms linear density fluctuation at the filtering scale. 
Various open symbols show computed power law slopes for 
all 25 256 3 models at four different redshifts, and respec¬ 
tive solid symbols show four 512 3 models at three highest 
redshifts. The solid line shows the fit in the form 




0.76 

ln(l + uf.) 0 - 9 ’ 


( 10 ) 


and dotted lines mark ±13% interval around the fit (which 
is consistent with being a 95% confidence interval). I con¬ 
clude that within the errors of the approximation, the slope 
of the column density distribution is a function of the rms 


linear density fluctuation at the filtering scale alone. This 
conclusion qu antita tively confirms the qualitative conclusion 
of Hui et al. (1997), who found that more nonlinear models 
have shallower column density distributions. 

The scatter is much smaller for the 512 3 models, as 
should be expected from Table However, since four 512 3 
do not cover the full space of cosmological parameters, they 
cannot be used to investigate whether the slope of the col¬ 
umn density distribution depends on other parameters in 
addition to of- I however use them to accurately determine 
the best 6 t. Fig. [I] thus demonstrates that the slope of the 
column density distribution of the Lyman-alpha forest f3 at 
jVh i, fit = 3 x 10 13 cm ~ 2 is the function of the density fluc¬ 
tuation on the filtering scale alone to within 13%. It is very 
likely, judging from how well points from 512 3 simulation 
fit the solid line, that /3 is a function of of to much better 
accuracy (at least 5%), but to verify this requires running 
several more 512 3 simulations. This work in now in progress. 

The scatter around the fit (|k]) may be caused by two 
different effects: by the error in the HPM-DPA calculation, 
and by dependence of the slope j3 on other cosmological pa¬ 
rameters, like the local slope of the power spectrum, the 
equation of state etc. Because the scatter is consistent with 
the error alone, current work is not able to address the sec¬ 
ond source of scatter, i.e. it is possible that the slope of the 
column density distribution does depend on other parame¬ 
ters of the model, but this dependence is weaker than 13% 
and therefore cannot be uncovered here. I would like to em¬ 
phasize that this conclusion does not necessarily contradict 
the conclusion of Hui et al. (1997), who found a strong de¬ 
pendence of the slope of the column density distribution on 
both the local slope of the density power spectrum and on 
t he eq uation of state slope 7 . The conclusion of Hui et al. 
( 1997 ) was based on using a quantity called 00 as a measure 
of the nonlinearity of the model. This quantity is ill-defined, 
because it measures the amplitude of linear density fluctu¬ 
ation on a scale which is model-dependent and lacks any 
physical explanation, while of represents to a high degree 
of accuracy (as have been shown in Gnedin & Hui 1997) the 
real total (i.e. summed over all scales) rms linear gas den¬ 
sity fluctuation. I also note here that since the Zel’dovich 
approximation used by Hui et al. (1997) is uncontrolled ap¬ 
proximation, its accuracy is not known. 

I also note here that the scatter around the fit gets 
larger for higher values of of, as could be expected, as higher 
of corresponds to more nonlinear models, for which the cos¬ 
mic variance is more significant. The crosses, which corre¬ 
spond to z = 2.31, have a larger scatter, and also somewhat 
biased toward higher values of f3. This is due to the fact that 
the range of column densities within which numerical con¬ 
vergence is achieved is not established at a = 2.31, as have 
been explained in the previous section. I show the simula¬ 
tion results at z = 2.31 here for illustration purposes only; 
they are not used to determine the fit ©■ 

Fig. ^ presents the slope (3 as a function of of, box, 
because it is of,box, rather than of, which is actually com¬ 
puted in a simulation. But now, given the dependence of 
the slope (3 on of, it is possible to correct for this bias, and, 
given the slope of the column density distribution derived 
from a simulation, /3sim, to compute the slope Pcor which 
would be derived in an idealistic simulation of an infinite 
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Peon = Psim + P(ctf) — P{ctf, box), 


( 11 ) 


and correct the column density distribution itself accord¬ 
ingly: 

IVh id 2 AT 


(1 + z)dz dNm 


N B 


(1 + z)dz dNm 


SIM 


N h i \ 


P(< 7 f)—P(. ct f, box) 


Nki. 


( 12 ) 


where Ahi,fit = 3 x 10 13 cm -2 . 

After having corrected for the bias due to the finite sim¬ 
ulation size, I can directly compare cosmological models to 
the observations. Let me first concentrate onz = 2.85, since 
for this value of redshift the HPM-DPA calculation is the 
most accurate. The obs erved slope of t he column density 
dist ribu tion is p = 0.5 (Hu et al. 1995; Kirkman & Tytler 
199? j) . h owever, the error of this number is not quite certain. 
Therefore, I will say that the model fits the observations if 
its slope p is between 0.4 and 0.6 (20% total [15% observa¬ 
tional plus 13% modeling] error); if the slope of the column 
density distribution for a model falls between 0.35 and 0.40 
or between 0.60 and 0.65 (30% total error), I will say that 
the model marginally fits the data; finally, if the slope P for 
a given model is above 0.65 or below 0.35, I will assume that 
the model fails. 

Table ^ presents the biased /3sim and corrected Peon 
values for the slope of the column density distribution for 
all COBE normalized models from Table [l|, as well as the 
value of p computed from equation (©■ The cosmological 
parameters for the models are included for reference. Finally, 
the last column indicates whether the model fits the data 
(+), marginally fits (±), or fails completely (—). 

Finally, using equation © , it is possible to place a con¬ 
straint on the value of the rms linear density fluctuation ap. 
The observational constraint 0.4 < P < 0.6 can be expressed 
in terms of of as 


1.6 <a F < 2.6. (13) 

This should be considered as roughly 2 a (95% confidence 
level) interval. With the 3 a error, this interval becomes 

1.5 <a F < 3.1. 


The median value, P = 0.5 corresponds to <jf = 2.0. 

Finally, as I have mentioned in the previous section, 
the filtering scale 1/Icf depends on reionization history and 
is a nontrivial function of time. For a purpose of a quick 
estimate, it would be more convenient to use a fixed scale 
like I//C 34 (eq. [(|). The following approximate expression 
accurate to about 5% can then be used to compute of at 
z = 2.85 from < 734 : 

a F (z = 2.85) = (7 3 4 [1.06 - 0 . 9(7 - 1.5) - 0.15(T 4 - 1)]. 

This expression depends on the parameters of the equation 
of state since (jf depends on the equation of state while 
(734 does not. I point out here that the above expression is 
derived by analyzing 25 models considered in this paper, and 
may not be valid for a model that falls outside the range of 
cosmological models considered here, or for a model whose 
equation of state evolves significantly d iffere ntly from the 
analytical expressions of Hui & Gnedin (|l998|). 



N«i (cm~ 2 ) 


Figure 5. Comparison between the observational data at z — 
2.85 and all 25 cosmological models, reduced to have the slope of 
the column density distribution of 0.5 (eq. [c | ) and normalized 
according to equation 11 4l (dotted lines). 


3.2 The amplitude of the column density 
distribution 

While the slope of the column density distribution can serve 
as a basis for discriminating against some of the cosmologi¬ 
cal models, the amplitude of the column density distribution 
depends on the assumed value for the ionizing intensity J 21 , 
which is currently highly uncertain. The best available mea¬ 
surements are based on the proximity effect, an d give values 
in the r ange J 21 = 0.2 — 1.5 at z = 2 — 4 ( GjaUongoet 
al. 1996 ; Lu et al. 1996| ; |Cooke, Espey, fe Carswell 1997j ). I 
therefore consider J 21 as a free parameter, and use the ob¬ 
servational data to fix the value for J 21 that best fits the 
data. 

Let me first consider z = 2.85. The normalization pro¬ 
cedure described in the previous subsection defines a value 
of J 21 for each cosmological model. After computing the re¬ 
spective J 21 values for all 25 models, I obtain the following 
approximate expression for J 21 as a function of cosmological 
parameters: 


J 2 i(z = 2.85) = 1CT 


fin /2 T„ 0 - 7 o 


(14) 


which is accurate to about 40%. The accuracy estimate re¬ 
sults form the assumed accuracy of 20 % for the amplitude 
of the column density distribution, which in turn translates 
into a 40% accuracy in J 21 oc Am for the power-law column 
density distribution with the index P « 0.5. 

In order to demonstrate graphically the level of accu¬ 
racy of this approximation, I plot all 25 models with the 
value of J 21 defined according to equation (|l4 ), “reduced” 
to have the same slope of the column density distribution, 

N Hl d 2 jV 


(1 + z)dz dAm 
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Table 4. Models vs observations 


Model 

n 0 

fib 

Qa 


h 

n 

^8 

PsiM 

PcOR 

P 

fit? 

SCDM.5A a 

1 

0.05 

0 

0 

0.5 

1 

1.2 

0.33 

0.32 

0.30 

— 

TCDM.2A 

1 

0.05 

0 

0 

0.5 

0.9 

0.83 

0.57 

0.51 

0.48 

+ 

LCDM.1A 

0.35 

0.05 

0.65 

0 

0.7 

1 

0.94 

0.59 

0.51 

0.51 

+ 

LCDM.1C 

0.35 

0.05 

0.65 

0 

0.7 

1 

0.94 

0.56 

0.52 

0.44 

+ 

LCDM.2A 

0.35 

0.05 

0.65 

0 

0.7 

0.96 

0.80 

0.63 

0.56 

0.55 

+ 

LCDM.3A 

0.35 

0.03 

0.65 

0 

0.7 

1 

1.04 

0.57 

0.51 

0.44 

+ 

LCDM.4A 

0.35 

0.03 

0.65 

0 

0.7 

0.95 

0.85 

0.65 

0.56 

0.58 

+ 

LCDM.5A 

0.4 

0.036 

0.6 

0 

0.65 

1 

1.02 

0.54 

0.48 

0.44 

+ 

CHDM.1A 

1 

0.05 

0 

0.1 

0.5 

1 

0.92 

0.52 

0.47 

0.44 

+ 

CHDM.2A 

1 

0.05 

0 

0.15 

0.5 

1 

0.86 

0.63 

0.54 

0.56 

+ 

CHDM.3A 

1 

0.05 

0 

0.2 

0.5 

1 

0.82 

0.72 

0.57 

0.73 

- 

CHDM.4A 

1 

0.05 

0 

0.1 

0.65 

1 

1.22 

0.39 

0.37 

0.31 

- 

CHDM.5A 

1 

0.05 

0 

0.2 

0.65 

1 

1.12 

0.56 

0.49 

0.46 

+ 

CHDM.6A 

1 

0.07 

0 

0.1 

0.6 

0.9 

0.73 

0.66 

0.56 

0.63 

± 


a A 512 3 simulation is used for this model. 


N B 


(1 + z)dz gWhi 


( AT \ 0.5-/3(t7 F ,BOx) 

\ 

-/Vh i,F] 


(15) 


In other words, if both equation (|H)]) and equation (0 were 
exact, the reduced column density distribution would be the 
same for all cosmological models. Thus, the difference be¬ 
tween different lines in Fig. |5| de monstrates the error in both 
equation ( |lo| ) and equation (|l4[). 

One can note that while the majority of models lie to¬ 
gether in a strip coinciding with the data points, there are 
six models (SCDM.2D, SCDM.2G, SCDM.2H, SCDM.2L, 
LCDM.1C and LCDM.2A) that lie above the rest. This dif¬ 
ference correlates with neither of the cosmological param¬ 
eters defining the model. More than that, since one of the 
outliers (SCDM.2D) is just a different random realization 
of the SCDM.2A model, the difference in the amplitude is 
real, and is due to the error induced by the HPM-DPA ap¬ 
proximation, most likely the finite size of a simulation box. 
One therefore may hope that the dispersion in Fig. ^ can be 
reduced by using larger simulation sizes (work currently in 
progress). 

Finally, I would like to point out here that fib and h 
dependence in equation ([Td]) agrees with the one obtained 
by Croft et al. ( 199f| ) when the dependence of To on fib (for 
a given cosm ological model) is taken into account ( Hui & 
Gnedin 1998). However, the advantage of equation (|1^) is 


that it is (approximately) valid for any cosmological model, 
i.e. there exist no additional dependences on physical param¬ 
eters that are not incorporated into equation (M), except at 
a level comparable to or below the error of the HPM-DPA 
calculation. 

I note here, that since both To and of are non-trivial 
functions of time, and it is unlikely that 2 = 2.85 represents 
a special moment in the history of the universe, equation 
( fL4| ) cannot hold at some other redshift with the coefficient 
lCr being replaced by some other appropriate number. In 
other words, while it is plausible that the expression ( floj ) 
for the slope of the column density distribution is accurate 
to high precision, equation (|IJ) is necessary approximate, 
and there has to be some other dependence in equation jl4| ) 
on cosmological parameters at the 40% or below level. 



Figure 6. Comparison between the models that fit the data at 
z = 2.85 and observations at four different redshifts. The same 
value of J 21 is adopted for each cosmological model at all four 
redshifts. Symbols show the observational data with solid symbols 
lying in the column density range used to fit the amplitude of the 
column density distribution. 


3.3 The evolution of the column density 
distribution 

Using the observational data at four different redshifts, it is 
possible to constraint the evolution of the column density 
distribution of the Lyman-alpha forest from 2 ~ 4 to 2 ~ 2. 

I show in Figure ^ the column density distributions for 
the seven models from Table |i] that fit the data at all four 
redshifts. For every model the value of J 21 is the same at 
each redshift (but different for different models). Since all 
models are consistent with the data at all four redshifts, 
the observations are consistent with the value of J 21 being 
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independent of redshift in the redshift interval z ~ 2 — 4. 
This normalization is also fully consistent with the slightly 
higher abundance of the Lyman-al pha f orest absorbers at 
z ~ 3 found by Kirkman & Tytler (1997) within the errors 
of the analysis presented in this paper. 


4 DISCUSSION 


I have demonstrated that the slope of the column density 
distribution of the Lyman-alpha forest is a function of the 
rms linear density fluctuation of the gas distribution (which 
is smooth on scales below the characteristic filtering scale), 
and is independent of other cosmological parameters within 
the accuracy of about 13% of the approximate method uti¬ 
lized (Hydro-PM approximation combined with the Density- 
Peak Ansatz). The amplitude of the column density distri¬ 
bution, expressed in terms of the ionizing intensity J21, is 
proportional to a combination of cosmological parameters, 


J21 oc 


rt /2 T°- V F ’ 


‘0 ± 4 

and is independent of other cosmological or physical param¬ 
eters within about 40% accuracy, at a given value of redshift. 

The observational data and the modeling technology 
are currently consistent with the value of J21 being constant 
from z 4 until z ~ 2. 

In order for a cosmological model to fit the observed 
slope of the column density distribution of the Lyman-alpha 
forest at z « 3, the rms linear density fluctuation at the fil¬ 
tering scale ~ /C34 (see eq. ||) in the model should be 
between 1.6 and 2.6 (or between 1.5 and 3.1 with a conser¬ 
vative estimate). 

Why is the rms linear density fluctuation the major fac¬ 
tor determining the slope of the column density distribution, 
when crp ~ 2 means an already nonlinear distribution? The 
answer is that of serves as a measure of the nonlinearity 
of the model; one could have adopted another measure, say, 
the total (nonlinear) rms density fluctuation in the model. 
However, since the total rms density flu ctuation is a function 
of the rms linear density fluctuation (Kofman et al. 1991) 
with, perhaps, a slight dependence on the local value of the 
slope of the density power spectrum, which is within the un¬ 
certainty of the simulations presented here, any measure of 
the nonlinearity of the model can be expressed as a function 


I compute the error in the amplitude of the linear density 
fluctuation at the cluster scale, 


08, CLUSTER 

The error in the slope of the column density distribution 
of the Lyman-alpha forest (assuming the observations give 
P = 0.5) is simply 



Table |E] lists all fourteen models together with the respec¬ 
tive errors in the density fluctuation amplitude at the clus¬ 
ter scale and in the slope of the column density distribution 
of the Lyman-alpha forest (the check mark after the error 
shows whether the model passes the test at a 95% confi¬ 
dence level or not). One can see that many models do well. 
Several of LCDM models pass all three tests. The fact that 
none of the CHDM models presented pass all the tests does 
not imply that this model can be ruled out: it seems plau¬ 
sible that a set of parameters can be found that will satisfy 
all three tests. The set of CHDM models used in this pa¬ 
per has been chosen “randomly”, without prior attempts to 
make these models fit the cluster abundance and the slope of 
the column density distribution tests. The fact that several 
of “randomly” chosen LCDM models perform well implies 
that there is more parameter space available for the LCDM 
models than for CHDM models, but it will be premature to 
claim that CHDM model fail altogether. 

It however seems possible to make a claim that the den¬ 
sity parameter in massive neutrinos, Q v , has to be lower than 
or about 0.15. The two models, CHDM.3A and CHDM.5A, 
represent two extreme cases of the = 0.2 models, and 
they clearly demonstrate that models with = 0.2 have 
too little power on scales of about 100 kpc as compared with 
lOMpc scales. 

I am very grateful to Romeel Dave for letting me use 
his AUTOVP automated Voigt profile fitting software. This 
work was supported by the UC Berkeley grant 1-443839- 
07427. Simulations were performed on the NCSA Power 
Challenge Array under the grant AST-960015N and on the 
NCSA Origin2000 mini super-computer under the grant 
AST-970006N. 


Of (TF- 

Finally, a few comments can be made about fitting a 
cosmological model to the observational data. If one ignores 
semi-philosophical questions about the Hubble constant, the 
age of the universe, and the value of flo, and restricts oneself 
to requiring that a given cosmological model satisfy three 
measurements of the underlying density power spectrum: 
the COBE normalization, the cluster abundance, and the 
slope of the column density distribution of the Lyman-alpha 
forest at z ~ 3, once can then ask a question how well the 
cosmological models listed in Table [l] would perform on those 
three tests. Let me concentrate only on the models that 
I the first test, i.e. only on the CORF, normalized models 


pass 


listed i: i Table Using the fit as ^miSXEB f° r the top-hat 
linear density fluctuation at the 8h _1 Mpc scale t hat sa tisfies 
the cluster constraint from Eke, Cole, & Frenk ( 1996| ), 

08,cluster = (0.50 ± 0.04)I2 o °- 53 + 013n o ^ 
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Table 5. CO BE normalised models 
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